Scaling of Non-Saturating MR and quantum oscillations in pristine and ion-implanted 

HOPG 
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A wide variety of resistive and field dependent behaviors have been previously observed in both 
doped and non-doped Highly Oriented Pyrolytic Graphite (HOPG). We find HOPG samples to vary 
significantly in their temperature dependent resistances, even between portions taken from the same 
sample, yet they exhibit consistent non-saturating magnetoresistance (MR). The scaling behavior 
of the MR is shown to be characteristic of a model based on the Hall effect in granular materials. 
In addition to the large, field-linear MR, all samples exhibit Shubnikov-de Haas (SdH) oscillations. 
Additional samples were doped via ion-implantation by boron and phosphorous, but show no signs 
of superconductivity nor any systematic change in their magnetoresistive behavior. Analysis of the 
SdH data gives a 2D carrier density in agreement with previous results, and a large mean-free path 
relative to crystallite size, even in samples with thin ion-implanted surface layers. 



Highly oriented pyrolytic graphite (HOPG) is a com- 
mercial product that finds uses in areas diverse as 
neutron monochrometers[lj and ultra-flat substrates for 
scanning probe microscopy. More recently, hints of 
magnetism and superconductivity have been reported in 
HOPG, the latter reportedly extending to 400 K[1|. Ion 
implantation of HOPG with phosphorous has been re- 
ported to induce superconductivity at low doses of low 
energy ions. Very large, non-saturating magnetoresis- 
tance and evidence for Shubnikov-de Haas (SdH) oscil- 
lations have also been found. [4i, ^] Large variations in 
transport and magnetic properties between samples, and 
even on portions of the same sample, are common. 
HOPG materials are typically characterized by their de- 
gree of mosaicity, by which is meant the angular spread 
of the c-axis orientation of component graphite crystal- 
lites. For example, ZYA quality is stated to have a mosaic 
spread of 0.4°, and is of the highest quality commercially 
available. Crystallite size (up to 10 fim for ZYA), inter- 
granular connectivity, and possible interstitial impurities 
are clearly candidates for the large variation in electrical 
and magnetic properties. 

In this paper, we focus on magnetoresistive (MR) ef- 
fects and attempt to separate intrinsic properties of the 
component graphite crystallites from the extrinsic effects 
associated with the composite nature of HOPG. We show 
that the microscopic properties, as revealed by the SdH 
oscillations, are independent of the bulk dc resistance. 
However, the very inhomogeneity of the material, com- 
posed of crystallites of different sizes and orientation, 
leads to important physical effects, with nonsaturating 
magnetoresistance [7] prominent among them. Similar 
observations of linear magnetoresistance have been re- 
ported on silver chalcogenides and were shown to persist 
to 60 T.fs^ An explanation was proposed by Parish and 
Littlewood, 9] who modeled the inhomogeneous conduc- 
tor as an array of coupled four-terminal circuit elements, 
as disks in which four currents can enter from adjacent 
elements. Each disk is characterized by its scalar resis- 
tivity p and Hall resistivity pfJ.H, where p is the mobility 
and H, the applied field. As the field increases, the Hall 



FIG. 1: Schematic structure of the typical sample of HOPG 
of ZYA quality used in experiments, with four contact leads 
attached to the surface, and the thin ion-implanted surface 
layer shown. The inset shows the orbits of electrons cross- 
ing the interfaces between crystallites perpendicular to the 
direction of the electric field. 



effect in each crystallite causes current to flow through 
terminals that are normal to the electric field, giving rise 
to the large linear MR that is observed. Allowing for 
an average mobility (p) and a standard deviation Ap, 
Parish and Littlewood determined the crossover between 
quadratic dependence on magnetic field, for fields be- 
low a characteristic value Hq, and linear behavior above. 
For Ap/ (p) small, as we might expect for HOPG, the 
crossover field is of order 1/ (p) and the overall scale of 
the MR, proportional to (p) . 

Samples for this study were obtained from ND-MDT 
Corp. and were of ZYA quality. Certain samples were 
subjected to ion implantation by boron at various dose 
levels and ion energies ranging from 10 to 40 keV in an 
effort to induce surface superconductivity within the ion- 
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FIG. 2: Temperature-dependent resistivity of a typical ZYA 
sample in zero field and in a field of 1 T on a logarithmic 
scale. The inset shows the zero- field data on a linear plot. 



FIG. 3: Contrasting zero-field resistance of a 0.5 mm thick 
sample of ZYA and that of a thin (0.05 mm) portion of the 
same sample 



implantation range. No evidence for superconductivity 
was found, and the MR behavior of irradiated samples, 
as will be seen, is identical to as-received material. With 
one exception, all samples were cut to approximately 2 
mm X 4mm, with four silver contacts placed across the 
width of the sample as seen in FigHJ The spacing of the 
inner pair (voltage contacts) was identical for all samples. 
Four-terminal resistivity measurements were carried out 
using a Quantum Design PPMS system over the tem- 
perature range 2K - 300 K and up to 7 T. The inset to 
Fig. 1 shows typical data on a ZYA sample which shows 
semiconducting- like behavior above « 100 K and metal- 
lic behavior below. As seen in the main portion of Figl^l 
the low temperature resistance increases by a factor of 
20 in a field of 1 T, and is semiconducting-like through 
the entire temperature range. 




As an example of the wide variation in the tempera- 
ture dependence of HOPG, Fig|3] shows the normalized 
resistivity a sample of thickness 0.5 mm along with that 
of a thin portion of the same sample having a thickness 
of 0.05 mm. While the thicker sample shows the same 
behavior as the sample in Figj^l the thinner sample is 
metallic over the entire temperature range. 

The nonsaturating MR is shown in Fig|3] for a sam- 
ple that had been ion implanted with boron at a level 
of 3 percent at an energy of 20 keV. At low tempera- 
tures, pronounced oscillations are observed which cause 
an apparent maximum in the resistivity near 15 K when 
plotted versus temperature at a fixed field of 3 T. Such 
downturns in resistance at low temperature have been at- 
tributed to the existence of superconducting regions, but 
are clearly a consequence of quantum oscillations. The 
magnetoresistance shows curvature at higher tempera- 
tures, but becomes increasingly linear as the temperature 



FIG. 4: The field dependence of the resistance at five tem- 
peratures for a sample of HOPG doped by boron ion- 
implantation. Note that Shubnikov-de Haas oscillations cause 
an apparent maximum in R(T) at 3 T, 



decreases. 

In FigEl we show the 5K MR of a ZYA sample to 7T. 
The linear field dependence is clear. The dashed line 
shows the data with the zero-field resistance subtracted. 
To the lowest fields measured, the data are not clearly 
quadratic in field 

As noted above. Parish and LittlewoodPl predict that 
the crossover from quadratic to linear field dependence 
occurs at Hq w 1/ (/z) when the variation of the mobility 
among crystallites is small. At the same time, the MR 
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FIG. 5: Logarithmic plot of the resistance at 5 K to 7 T. The 
mgnetoresistance is linear in field (red line) from 0.7 T to 7T, 
with no evidence for saturation. The dashed curve shows the 
same data with the zero-field value subtracted. The data are 
not purely quadratic in field down to 0.01 T. 
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FIG. 6: Left figure: resistance curves at four temperatures 
scaled by a characteristic field that follows the inverse of 
the mobility, decreasing with decreasing temperature, with 
inset showing the dependence of i^o on temperature. Right 
figure: the same data before scaling. 



itself is proposed to scale with Hq. We test this prediction 
in Figini at temperatures where the MR oscillations are 
less evident. is chosen for each curve at the point 
where it deviates from linearity by a set percentage. Note 
that the crossover field increases with temperature, 
reflecting the decrease in mobility as expected for phonon 
scattering in each graphite crystallite. The cross-over 
field extrapolates to = 0.08T at low temperatures. 

We turn now to the low temperature magnetoresis- 
tance oscillations. Early work by Williamson et al.p^ 
on the de Haas-van Alphen oscillations of pyrotylic and 
single crystal graphite showed their Fermi surface in 
detail. The majority electron and hole carriers were 
found in approximately ellipsoidal pockets along HKH 
in the Brillouin zone. There are six such ellipsoids with 
one third of their volume within the unit cell. Previ- 
ous work on Shubnikov-de Haas oscillations in graphite 
shows that the majority electron ellipsoids dominate in 
the resistivity. [Tl| 

Despite wide variations in temperature dependence of 
various samples, the MR is surprisingly similar. The 
main panel of FiglT] shows the numerical derivative 
dR{H) /dH of the data in FiglS] plotted vs the inverse 
field. This method of separating the oscillatory behav- 
ior from the linear MR was reported by Arnt, et a\.fd\ 
The periodicity in identifies these as Shubnikov- 

de Haas oscillations within each graphite crystallite. The 
inset of Figl7]shows similar curves for both doped and un- 
doped samples whose zero field resistance can vary by a 
factor of two. We note that the samples have the same 
thickness and contact spacing, so that the data actu- 
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FIG. 7: The data for HOPG ZYA (points), shown with the 
fitting curve (line). The inset shows the Shubnikov-de Haas 
oscillations for the ZYA as well as two ion-implantation dop- 
ing levels of 3 and 4 percent. 



ally reflects nearly identical magnetoresisistivity, includ- 
ing MR oscillations. 

Fitting the curves in FigjT] to a version of the SdH 
oscillations with a constant relaxation timejl2j, 



dR 

d^oH 



27r/3 



(1) 



yields the linear MR coefficient Rf = 0.084r2/T, oscilla- 
tion amplitude A = 0.18f2/T, damping factor a — 2.7T; 
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period /3 = 4.6T and phase factor (j> = — 0.127r. The 
standard expressions [l^ [l3j for Shubnikov-de Haas os- 
cillations give 



where Sp is the cross sectional area of the Fermi sur- 
face. We obtain kp = {Sf/ttY''^ = 1.2 x 10^m-\ which 
leads us to a 2D carrier density of 

n2i, = ^ = 2 X IQi^m-^ (3) 

The Fermi surface for holes can be seen weakly in de 
Haas- van Alphen oscillations but not SdH.fTo^ 

Studies of microcrystalline graphite [5] find fc^? « 6 x 
lO'^m-i with a 2D carrier density ?s 6 x 10^** m ^. That 
value agrees with our result to within a factor of 2. We 
note that each Landau Level contains a 2D density of 
states of -ff/2$o = (-ff/4T) x 10^^ m^^^ where $0 is the 
flux quantum. At the first peak in Figl?! where H — A 
T, only the first and second Landau levels are filled, each 
with 10^^ m~^. As the field is increased through 4 T, 
the second Landau level begins to empty and there are 
no further oscillations. 

Treating the majority electron pocket as an ellipsoid 
of radius kp with a height ^ that of the Brillouin zone, 
we arrive at a 3D carrier density of 2.3 x 10^'*m~^. This 
is in reasonable agreement with the value given by band 
theory of 3.3 x lO^^m-^.fl^l 

The damping factor provides a direct measure of the 
carrier mobility through a — n/fj, 01 jj, — 1.2x10^ 
cm^/Vs. This corresponds to a mean-free path « 0.1 /im, 
indicating that the crystallites contributing to the S-dH 
oscillations and the large linear MR are substantially de- 
fect free. Expressed in field terms, the mobility is fi = 
1.2 T~^. It is also important to remember that this mo- 
bility is an average over individual crystallites, and not 
that of the bulk. If, as predicted, the cross over between 



quadratic and linear MR regions is inversely proportional 
to fi, we would expect Hq = 0.8 T. In fact, the scaling 
fields in Figl6] extrapolate to a low-temperature value of 
0.08 T. In the Parish-Littlewood model, the crossover can 
be dominated by the spread in mobility, but this would 
require that A/i ~12 T^^, which does not seem reason- 
able. Perhaps further refinement of that model will pro- 
vide the needed proportionality factor (J-Hq ~ 0.13. 

In conclusion, while there is a significant variability in 
the temperature dependence of the resistance of samples 
without an applied field, even between portions taken 
from the same sample, and particularly for ion-implanted 
samples, the large non-saturating magnetoresistance is 
consistent throughout. This MR is characteristic of the 
Parrish-Littlewood model where the currents flow in a 
nonhomogeneous system of graphite crystallite disks per- 
pendicular to the electric field in the presence of an ap- 
plied magnetic field as shown in Fig[TJ The numeri- 
cal simulation in the model predicts that AR « R{H — 
0)H/Ho, when the cross over is dominated by the mobil- 
ity From Eq. ^ and Fig [5] we have i/o = 0.08 T and 
R{H = 0) = 4 mfJ which predicts AR = 0.35 fl at 7 
T, which compares well with the data in FiglS] We note 
that saturation is not predicted to occur in this model 
until fiH Z 20; our value at 7 T is only ^iH = 8.4. The 
SdH oscillations are in agreement with previous results 
for the majority-electron Fermi surface and give a major- 
ity electron density consistent with band theory. The fact 
that only the majority electrons contribute to the SdH 
oscillations, allows us to calculate the mean free path, 
showing the crystallites comprising the HOPG sample to 
be relatively defect free. It means that the ion-implanted 
region that has many defects due to the impact of ions 
does not contribute to the SdH oscillations, implying that 
the current distribution does spread far below the ion- 
implanted surface region so that the average mobility is 
not significantly changed by ion-implantation. No evi- 
dence of ion-implantation induced superconductivity has 
been detected for the ion energies and does used. 
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